Two-dimensional wave patterns of spreading depolarization: 
retracting, re-entrant, and stationary waves 
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Abstract 

We present spatio-temporal characteristics of spreading depolarizations (SD) in two experimental systems: retracting SD wave 
segments observed with intrinsic optical signals in chicken retina, and spontaneously occurring re-entrant SD waves that repeatedly 
spread across gyrencephalic feline cortex observed by laser speckle fiowmetry. A mathematical framework of reaction-diffusion 
i systems with augmented transmission capabilities is developed to explain the emergence and transitions between these patterns. 
Our prediction is that the observed patterns are reaction-diffusion patterns controlled and modulated by weak nonlocal coupling. 
The described spatio-temporal characteristics of SD are of important clinical relevance under conditions of migraine and stroke. 
In stroke, the emergence of re-entrant SD waves is believed to worsen outcome. In migraine, retracting SD wave segments cause 
neurological symptoms and transitions to stationary SD wave patterns may cause persistent symptoms without evidence from 
noninvasive imaging of infarction. 

Key words: pattern formation, reaction-diffusion, spiral waves, migraine, stroke 
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1. Introduction 

The World Health Organization lists migraines among the 
four most disabling chronic medical disorders Oh- Stroke is 
currently the third leading cause of death in developed coun- 
tries. There is clinical and experimental evidence that in both 
conditions spreading depolarization (SD) waves of non-or-all 
type emerge. SD waves are characterized by massive redis- 
tribution of ions across cell membranes peaking after several 
seconds in a nearly complete neuronal depolarization, followed 
by a much slower recovery process taking up to minutes dur- 
ing which ion gradients are re-established towards their physi- 
ological values. This maximal ionic perturbation clearly distin- 
guishes SD from all other brain states such as epileptic seizure 
activity, functional activation or the physiological resting state. 
The sequence of ionic perturbation and its recovery spreads at 
a pace of about 3 mm/min over cortical regions. 

In migraine, this phenomenon — originally termed spread- 
ing depression of activity in the cerebral cortex J2l, causes neu- 
rological symptoms starting about 30 min before the headache 
phase |H E[]. It is also currently debated whether migraine 
headache is caused by substances released in the course of SD 



jfj. In stroke, so-called periinfarct depolarizations J3l, a vari- 
ant of SD, occur and spread in the boundary zones surrounding 
a cerebral infarct. A greater number of SD waves is believed to 
worsen outcome. SD has been demonstrated in patients suffer- 
ing from traumatic brain injury, ischemic stroke, subarachnoid 
haemorrhage and intracerebral haematoma 0,0, HE. Ill- 
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The difference between SD in migraine and stroke is that 
SD in the latter occurs when the cortex is deprived of ade- 
quate perfusion, i. e., under ischemia, while SD in migraine oc- 
curs under normoxic conditions, but otherwise both processes 
share a common underlying neuronal mechanism Il2ll . and they 
should not be seen as two different processes. For this reason 
they are referred to as spreading depolarizations 

03]. This 

name reflects the common mechanism rather than the effect (de- 
pression of activity) or location (periphery of infarct) of these 
waves as the other terms do. To focus on common macroscopic 
features is in particular advisable when two-dimensional spatio- 
temporal characteristics are investigated. Moreover, it empha- 
sises the depolarization wave as one of the common aspects 
within the complex bidirectional relation between migraine and 
stroke 031 [l4ll . For example, SD activity may play an im- 
portant role in complications of migraine, including migrain- 
ous infarction, persisting migraine aura without infarction, and 
ischaemia-induced migraine (symptomatic migraine). 

In this study, our aim is to understand the emergence of 
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pathological spatio-temporal SD states in the cortex and pos- 
sible transitions between them in terms of nonlinear dynam- 
ics and pattern formation. This approach complements more 
traditional clinical research methodologies. The paper is or- 
ganized as follows. In Sec. [2] we introduce the experimental 
setup, followed by a description of the observed patterns in 
Sec. [3] In Sec. |H we propose a mathematical framework to 
model reaction-diffusion patterns that are controlled and modu- 
lated by weak nonlocal transmission, and in Sec. [5] we illustrate 
the applicability of this concept by special cases derived from 
this framework. We end with a discussion in Sec. [6] 

2. Methods and Materials 

Monitoring the spatio-temporal characteristics of SD in two 
spatial dimensions is crucial for studying this phenomenon, in 
particular from a clinical point of view. We present data from 
two experimental systems of complementary complexity each 
with a specific optical method. Optical methods have much bet- 
ter spatial resolution than electrical measurements. They pro- 
vide indirect evidence of SD. The information about the spatio- 
temporal development of SD is very precise since correlation 
between the optical signals and electrophysiological changes of 
SD are well studied IU5l fl6[ [l7 1 . 



2.1. Chicken retina as SD model 

The chicken retina is a typical in vitro system used to study 
the formation of pathological activity pattern of SD by intrinsic 
optical signals lfl8l[l9ll . Retinal SD can be easily observed by 
intrinsic optical signals, because it evokes an optical signal with 
an amplitude being several orders of magnitude higher than that 
observed in cortex. In fact, retinal SD can be observed by the 
naked eye as a milky front. Moreover, the chicken retina is a 
very homogeneous neural tissue because it lacks blood vessels. 
The neural aspects of pattern formation can be thus studied in 
retinal SD in isolation, that is, without the neurovascular cou- 
pling, with the obvious disadvantage that retinal SD provides 
no information on this coupling. Last but not least, the retina 
has, similar to the cortex, a layered structure with several lay- 
ers of neurons interconnected by synapses and there is a solid 
knowledge of retinal network functioning I20I1 . 

The chickens (4-7 days) were decapitated and the eyes were 
carefully extracted from the orbit. The eye was purged from 
adherent muscles and connective tissue and afterwards the eye 
was equatorially divided. The vitreous body was gently re- 
moved, and then the eyecup including the retina was transferred 
into a small petri-dish filled with standard Ringer (R) solution 
(Fig 03 (a)). The standard Ringer solution has the following 
composition (in raM): 100 NaCl, 6 KCl, 1 CaCh, 1 MgS0 4 , 
30 NaHC0 3 , 1 NaH 2 P0 4 , 30 glucose. Throughout the ex- 
periment, the submerged retina was perfused with fresh Ringer 
solution. For more details see Ref. [18]. 

Following this preparation protocol, the healthy chicken retina 
has a dark appearance. Since blood vessels do not course across 
the retinal surface, the only landmark, situated in the infero- 
temporal field, is the pecten oculi (FigQ](b)), a highly vascular- 
ized organ with corrugated structure that projects from the optic 



disc into the vitreous body of the avian eye. The pecten is be- 
lieved to nourish the retina. Pectineal blood capillaries deliver 
nutrients to the retina, therefore, there are no retinal blood ves- 
sels. This fact, together with the intrinsic optical signal, makes 
the avian retina an ideal experimental system to study SD. 

All measurements were obtained with a CCD video camera 
and stored with a DVD-recorder. The retina was pricked with a 
fine Tungsten needle to elicit a retinal SD wave. After measur- 
ing a control wave, the standard Ringer solution was replaced 
with a Ringer solution where drugs had been added. After a 
time of 30 minutes a new wave was elicited. The experiments 
in chicken retina were approved by the local ethics commit- 
tee (Stuttgart and Magdeburg) and the Regierungsprasident of 
Stuttgart and Magdeburg. 

2.2. Stroke model 

The choice of the cat as a species for studying SD has been 
governed by disappointing experience with the use of rodent 
models of stroke: none of the promising trials of neuroprotec- 
tive agents in rodents, especially those that block SD, has been 
reproduced in clinical trials of the same agent. The cat cortex 
is closer than rodent cortex to that of humans in two impor- 
tant respects. First, it is gyrencephalic (folded), whereas the 
rodent cortex is lissencephalic (smooth). The data obtained es- 
tablish beyond doubt that SD spreads via sulci (clefts between 
gyral convexities), and never "jumps" between gyral convexi- 
ties; the information has allowed reliable interpretation of data 
from monitoring the human brain, and indeed provided a justi- 
fication for the initial work (with invasive monitoring methods) 
seeking evidence of depolarisations in the injured human brain. 
Secondly, the ratio of astrocytes to neurones in the cat brain is 
much closer to that of the human brain than it is in the rodent 
brain. This is important since the small energy metabolite pool 
in the brain, glycogen, is almost entirely located in astrocytes, 
again making a model of occlusive stroke in cats a more reli- 
able analogue of human stroke than one in rodents. A further 
advantage of the cat brain is that the anatomical arrangement of 
its blood supply closely resembles that of the human brain. 

In cat cortex, SD can be visualized by illuminating the cor- 
tex with laser light and imaging the resulting speckle pattern. 
The laser speckle contrast imaging can be used as a proxy for 
SD, because cerebral blood flow changes are almost invariably 
coupled to SD and so these changes can be studied as sur- 
rogate markers to image and map the occurrence and spread 
of SD in real time I! 611 . The imaging observation field was 
chosen so as to achieve maximal spatial resolution in the area 
around the periphery of the initial stimulation site. Fig. [3(a) 
illustrates the anatomy of the cat cerebral cortex, where the 
gyri are arranged in concentric arcs, in the order ectosylvian 
(EG), suprasylvian (SG), marginal (MG) gyri. The experiments 
in cats were approved by the local ethics committee and the 
Regierungsprasident of Cologne. 

3. Experimental results 

In the following two subsections, we present the observed 
SD patterns in two experimental systems. We then develop in 
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Sec.|4]the conceptual framework for modeling SD on a macro- 
scopic level that allows us to understand these experimentally 
observed spatio-temporal patterns of SD and possible transi- 
tions between them in terms of nonlinear dynamics. In Sec. [5] 
special cases are derived from this framework. 



3.1. Retracting retinal SD waves segments 

An SD wave segment is a wave with a broken wave front 
and thus with two freely moving open ends. The open ends of 
the wave front allow more complex spatio-temporal dynamics, 
in contrast to a wave with a closed wave front or a wave where 
the open ends are attached to a non-excitable boundary. Moni- 
toring the dynamics of the open ends can be used to probe the 
tissue excitability, as we will first show experimentally and then 
develop the mathematical framework in the next section. 

Travelling SD wave segments have, to our knowledge, nei- 
ther been observed in cortex nor in any other grey matter brain 
region of any species. However, indirect evidence, presented in 
a recent study combining migraine patients' reports and high- 
field functional magnetic resonance imaging, suggests that re- 
tracting SD wave segments cause migraine aura symptoms 121 1. 
We show here that these retracting waves can be observed in 
retinal SD by controlling excitability in a regime weakly sus- 
ceptible to SD. This provides the first direct evidence that the 
recently predicted retracting waves exist in SD and that their 
emergence depends critically on the tissue excitability. 

An SD wave segment can be created from an initially closed 
circular SD wave front. The circular SD wave is created by lo- 
cal superthreshold stimulation. One way to break the SD front 
is by manually inducing a transient excitation block in the tis- 
sue in some confined location in front of the approaching wave. 
Performing this procedure, the two open ends of the wave seg- 
ment usually curl in to form a double spiral, whose wave front 
then repeatedly re-enters the same retinal tissue. In chicken 
retina, a spiral-shaped SD can rotate up to 20 times before its 
center collides with the tissue boundary 12211 . Another way to 
obtain wave segments is to reduce retinal susceptibility to SD to 
a degree at which the wave front spontaneously breaks. In this 
case, the newly emerged open ends retract, as we show, and the 
wave vanishes invading the tissue only once. 

Firstly, we describe the manual creation of a spiral-shaped 
SD wave to contrast this procedure and the resulting spatio- 
temporal development with retracting SD wave segments. Reti- 
nal SD is elicited mechanically by gently touching the tissue 
with a sharp needle. In Fig.Q](c), the stimulation needle (S) is 
seen by its dark shadow. After about 15 s, when the extension 
of the circular SD wave front has grown in diameter to about 
0.75 mm, we administer locally through a pipette a 0.5 ml drop 



of Ringer solution containing a tenfold raised Mg 2+ concentra- 
tion ([Mg ]r = 10 mM). Pipettes having an inner tip diameter 
of about 0.5 mm are best suited. They must be placed directly 
over the surface of the submerged retina without touching it and 
less than 1 mm in front of the approaching SD wave (Fig.Q](e)). 
For visual guidance, the contours of the pipettes, which is out 
of focus in Fig.Q~](e), is indicated by dashed lines and its target 
position marked by an asterisk. The propagation of the circular 
SD wave is blocked around the pipette's target position and as 




Figure 1 : Creation of a retinal SD wave segment with freely revolving open 
ends: (a) scheme of experimental set up with eyecup preparation placed in a 
petri-dish, perfused with Ringer (R) solution and monitored by a CCD camera; 
(b) top view image from eyecup; (c) mechanical stimulation with sharp needle 
S; (d) circular SD wave, (e) local administration of Mg 2+ via pipette indicated 
by dashed lines at location marked by an asterisk; (f) retinal SD wave segment. 



a consequence its wave front breaks open. This wave segment 
develops into a spiral-shaped wave (Fig.[T|(f)), if the submerged 
retina is kept in standard Ringer solution and one end attaches 
to the tissue boundary of the retina. 

We now adjusted retinal susceptibility to SD by adding co- 
caine in a 1 - 3 mM concentration range in 1 mM concentration 
steps. At 1 mM cocaine concentration, SD waves propagated 
normally after being stimulated mechanically. At 2 mM cocaine 
concentration, we observed a spontaneous excitation block. In 
Fig. |2] (c), this happens at a location where the retina appear- 
ance was less dark indicating a compromised tissue state. The 
tissue at this location is, however, still viable. After the cocaine 
is washed out, SD waves propagate normally again in the en- 
tire retina. The open ends of the SD wave segment retracts in 
the further course (Fig. [2] (d)-(f)). In Fig. [2] (b)-(f), we man- 
ually and individually enhanced the brightness and contrast in 
rectangular regions of interest to better visualize the SD wave. 
Furthermore, in figures (Fig. [2] (g)-(l)) the spatio-temporal de- 
velopment is schematically shown in a drawing. 

The retinal susceptibility to SD was at 2mM cocaine con- 
centration already so low, that in most cases the initiated wave 
collapsed after spreading only a short distance. The term col- 
lapse refers to an unstable, and therefore transient, SD wave 
profile that fades away although the front does not break. This 
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Figure 2: Six snapshots taken at 40 s time interval from chicken retina, (a)-(f) 
original recordings from CCD camera with enhanced brightness and contrast 
in rectangular regions of interest. White arrow in (c) marks the breaking of 
the retinal SD wave, and subsequent pair of white arrows in (d)-(f) mark the 
location of the open ends of retinal SD wave, (g)-(l) Schematic illustration of 
this sequence. Retinal SD wave shown as black line, tissue in a compromised 
but viable state is marked gray, and pecten is shown filled with stripe pattern. 



process was predicted for SD to occur extremel y cl ose to the 
susceptibility range of retracting waves segments [23 ]. At 3mM 
cocaine concentration retinal SD waves could not be elicited 
anymore. Retracting retinal SD waves have also been observed 
by reducing the potassium concentration in Ringer from [K + ]r - 
6 mM to values between 1-2 mM. For concentrations below 
[K + ]r = 3 mM, the intrinsic signal of the retinal SD wave pro- 
file changes in a characteristic way II24I1 . and the condition of 
the retina was hardly stable for more than 30 min. Therefore no 
control measurement perfusing the retina with standard Ringer 
could be done. 



3.2. Re-entrant cortical SD waves 

A single SD wave that repeatedly invades the same tissue 
forms a re-entrant pattern. A re-entrant pattern is easy to under- 
stand as a wave propagating around an obstacle. This obstacle 
can be a functional or anatomical block. In this subsection, we 
present re-entrant SD patterns revolving around an anatomical 
block that develops after microinjections of high concentration 
of potassium have been locally administered to the cortical sur- 
face. Firstly, we describe further in the next paragraph func- 
tional blocks in SD and re-entrant spiral patterns, whose cre- 




Figure 3: (a) Diagram of the lateral aspect of the right hemisphere of the cat 
brain illustrating typical gyral anatomy of a small ischaemic area (gray) within 
the area exposed for imaging (rectangle), (b) Following KCL injection (at t=0), 
and an immediate transient radially spreading wave from the injection site, a 
new single wave of hyperemia appears spontaneously and cycles anticlockwise 
3.5 times around an injured area of some 3 - 5mm diameter centred on the 
original injection site, and commencing medial to it (above as viewed in this 
image field). Two snapshots were taken at a 3.5 min interval. The reader is 
advised to run the video to see the patterns. 



ation protocol was introduced in the previous subsection (Fig.Q] 
(f)). 

A functional block refers to a block that transiently forms as 
an intrinsic property of the 2D wave dynamics. The core of the 
rotating spiral center in retinal SD is an example of a functional 
block. Since re-entrant SD waves with a functional block in the 
cortex have, to the best of our knowledge, not been reported, 
we refer here to the observed characteristics in the retinal ex- 
periments. The spiral core can be a circular area under certain 
conditions, but it is in general a more complex structure. For ex- 
ample, the core does not need to cover a spatially extended two- 
dimensional area, instead it can be an one-dimensional curved 
line (inset (a) Fig. [5j, which itself performs a rotation, as is the 
case for retinal SD I22I1 . The detailed dynamics of the core is 
a measure of susceptibility to SD and provides important infor- 
mation about macroscopic model parameters (see Sec. 14.6b . 

In contrast to a functional block, an anatomical block refers 
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to non-excitable tissue, possibly an infarct region or an artifi- 
cially created lesion. SD patterns rotating around an artificially 
created lesion have also been called reverberation of SD waves 
1 25]. Reverberation of SD was demonstrated in rat cerebral cor- 



tex and has been known for a long time. Recently, data from 
monitoring the human brain has renewed interest in re-entrant 
patterns II ill , in particular in experimental models of stroke to 
investigate the cause of infarct expansion, because more fre- 
quent SD waves are believed to worsen stroke outcome and 
therefore the emergence of re-entrant SD waves is of consid- 
erable clinical interest. 

The creation of a re-entrant wave demands, in general, spe- 
cial initial conditions. SD is induced with microinjections into 
the cortex of 1 M potassium chloride (KCl). As described in 
Sec. 12.21 laser speckle contrast imaging flowmetry is used to 
detect waves of regional blood flow alteration coupled to SD 
in a cortical field of view covering three gyri called marginal, 
suprasylvian, and ectosylvian gyrus (MG, SG, and EG) respec- 
tively (Fig.[3](a)). A radial pattern is regularly observed for the 
initial wave after an injection, starting from the point of KCl 
injection and propagating outwards, often appearing as a con- 
centric wave. This is demonstrated with a supplementary video 
(Video 1, filename dahlemEtAl Video. mov). 

Following the radial pattern, an area of markedly reduced 
perfusion becomes established surrounding the injection point. 
We infer that an ischemic region develops at this site, most 
likely due to the marked vasoconstriction associated with high 
levels of K + concentration in the perivascular fluid. It is known 
that the perivascular administration of K + concentrations at lev- 
els of 40 mM decreases pial arteries calibre by 47.2% and vein 
calibre by 13.5% l26ll . This significant reduction changes the 
conditions at the injection point and in its border zone. At least 
this assumption is consistent with the changes following the ra- 
dial pattern and being observed in the surrounding of this de- 
veloping ischemic focus. 

During observation for 40 minutes following a KCl injec- 
tion, we saw new waves of hyperemia that developed sponta- 
neously at the edge of the lesion, without any further KCl in- 
jection. Surprisingly, these secondary waves did not propagate 
radially but instead each revolves several times the ischemic fo- 
cus. This is demonstrated in the supplementary video (Video 1, 
filename dahlemEtAlVideo.mov). Similar perilesion re-entrant 
patterns have been seen around experimental ischaemic lesions 
in the cerebral cortex after occlusion |l7l| . Furthermore, note 
that re-entrant patterns result in consistent periodicity of depo- 
larizations when monitored at a single cortical point, and such 
peri odic events have also been seen in the injured human brain 

4. Augmented reaction-diffusion model 

In this section, we develop the conceptual framework for 
modeling SD on a macroscopic level. This framework consists 
of a continuum model of the cortex in which the communication 
between neurons is separated in two major schemes: (i) local 
coupling, and (ii) nonlocal coupling in space and time. In the 
cortex, these schemes represent (i) diffusion, and (ii) network 



communication and neurovascular feedback. With this frame- 
work general information in gained on how reaction-diffusion 
patterns of phathological states in the cortex are controlled and 
modulated by weak nonlocal transmission schemes. We ap- 
ply this framework in Sec.[5]to interpret the experimentally ob- 
served patterns, in particular, their emergence and transitions. 

4.1. Model equation 

The level of detail and complexity of our framework is de- 
termined by both spatial extension and propagation speed of the 
emerging patterns. SD waves extend in the cortex over several 
centimeters and spread with a remarkably slow speed of several 
millimeters per minute. The extension of SD patterns suggests 
that its mathematical description should be in terms of large- 
scale activity in neural populations and describe how this ac- 
tivity affects the cortical ionic equilibrium (homeostasis). The 
slow speed of SD (about 3mmmin = 50 fj.ms~ l ) indicates 
that diffusion in cortical tissue plays a major role. Therefore, 
the central parts of a macroscopic continuum model of SD are 
reaction rate terms and diffusion terms. These terms are supple- 
mented by additional terms reflecting augmented transmission 
capabilities. Within such a framework, the model equation is 



d t u 

e~ l d t v 



f(u,v) 
g{u, v) 



+ D 



+ H 



(1) 



where d, is a short hand notation for the derivative with respect 
to time. The terms on the right hand side represent (i) reaction 
rates f(u, v) and g(u, v) of a simplified activator(«)-inhibitor(v) 
model of cortical homeostasis (Sec. 14.21 ). (ii) diffusion of chem- 
ical signals and other local transmission schemes summarized 
under the term volume transmission in the brain and described 
by the diffusion operator D (Sec. 14.31 ). and (iii) the augmented 
transmission capabilities described by an operator H (Sec. 14.41 ). 
The model equation ([]]) and also its derived special cases (Sec. 
[5]) are nondimensionalized equations. In such systems, natural 
units of a system are used in a systematic manner so that the in- 
put and output variables of the rate functions f(u, v) and g(u, v) 
are in the order of unity. Therefore, beside many other advan- 
tages, nondimensionalization suggests explicitly a parameter, 
namely e that should be used for analyzing this pattern form- 
ing system. The parameter s indicates how much slower the 
production rate of the inhibitor v is compared to the one of the 
activator u. 

4.2. Activator-inhibitor dynamics 

Local dynamics of activator-inhibitor type are common in 
pattern forming systems, and — if they are coupled by diffusion- 
have been shown to gov ern several aspects of SD wave prop- 
agation J23, IH, |29l M, 31]. The reaction rates f(u,v) and 
g{u, v) in Eq. (Q]) represent such simplified local dynamics in the 
cortical ion homeostases (see Fig. HI leading to non-or-all type 
behavior (excitability). The rate function f(u, v) describes the 
change in a proxy agent named activator u by lumping together 
nonlinear reactions containing a positive feedback loop in the 
SD mechanism, such as A^a + -inward currents (dashed arrows 
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Figure 4: Flow chart of feedback loops and transmission pathways in SD. 

in Fig.|H and extracellular potassium concentration [K + ~\ e (dot- 
ted arrows). Likewise, the rate function g(u, v) of the inhibitor 
v represents recovery processes, such as effective regulation of 
[K + ] e and pump currents by the neuron's electrogenic Na + -K + 
ion pump (solid arrows with circular heads). 

An activator-inhibitor system exhibiting excitability requires 
a certain configuration of trajectories that follow the vector field 
give n by the functions f{u, v) and g{u, v) in the phase space 
1 321, L33Q - This configuration usually results from the parameter 
vicinity of an oscillatory regime whose large amplitude limit 
cycle is suddenly destroyed. A rest state (fixed point) is usu- 
ally the only attractor in the excitable regime. There exists 
also a threshold close to this rest state, e.g., a separatrix or a 
sharp trajectory, that must be crossed by a stimulation to result 
in phase space excursion. Activator-inhibitor systems are the 
most generic systems to exhibit such behavior. A particular ex- 
ample is derived in Sec. l4.5l and the validity and applicability of 
generic dynamics to model SD is discussed in Sec. [6] 

4.3. Volume transmission 

Communication between neurons is usually classified in two 
schemes: wiring and volume transmission [34]. Volume trans- 
mission is mainly characterized by diffusion of chemical sig- 
nals in the extracellular space. That diffusion is the main cause 
of the spread in SD is indicated by the propagation speed of 
about 50 11m s~ l . Waves caused by nonlinear reactions coupled 
by diffusion of monovalent ions travel at speeds in the order 
of lO/zms -1 assuming moderate reaction rates in the nonlinear 



dynamics Ir3lll . Signals transmitted by cortical wiring, in par- 
ticular action potential along axons, travel at speeds in the order 
of meters per second, i. e., five orders of magnitude faster. 

The diffusion operator in Eq. (|TJ and takes the general form 



D 



1 
A 



(2) 



with V 2 being the Laplace operator. We do not consider cross- 
diffusion, therefore off-diagonal entries vanish. The foam-like 
geometry of the interstitial space renders diffusion in the brain 
different from aqueous solution. In particular, apparent dif- 
fusion coefficients should be considered in a dimensionalized 
model Ir35tl . In our nondimensionalized framework, diffusion 
coefficients of activator and inhibitor, D u and D v , respectively, 
do not appear explicitly. Only the ratio A = of the diffusion 
coefficients enters Eq. (f2]i as a parameter (note that s appears 
also in this quotient). In SD, inhibitory reactions are thought to 
relate to ion pump activity (Fig. 0), in which case the inhibitor 
would not be subject to diffusion and therefore we choose A — 
in this study. 

4.4. Augmented transmission 

Similar to the diffusion operator in Eq. ©, the operator H 
describes the augmented transmission capabilities as 



I A "" A A m I 



(3) 



F is a operator specifying the augmented transmission (in anal- 
ogy to the Laplace operator V 2 for volume transmission), ex- 
amples are given in Sec. The four elements of the matrix Ay 
specify the relative coupling strength among different schemes 
of augmented transmission (in analogy to A in Eq.[2ji so that the 
operator AjjF represents three individual steps of augmented 
transmission, namely (I) selecting a species j whose transmis- 
sion capability is augmented, (II) creating the nonlocal driving 
force from this species, and (III) feeding this driving force back, 
multiplied by the coupling strength K, into the dynamical vari- 
able of species ;, with i, j e {u, v}. 

The parameter K describes the coupling strength relative to 
the diffusion coefficient, which is normalized to unity. There- 
fore, the Ajj are also normalized to unity, i. e., 



max({|A J, \AJ, \A VU \, \A VY \}) = 1 
and we limit our framework to cases where 
<1 



(4) 



(5) 



holds. With the constraints in Eqs. (IUi-© diffusion becomes 
the dominant coupling. We pose these constraints, because, 
within this framework, augmented transmission capabilities must 
not become an essential part of the SD wave propagation mech- 
anism. Rather we assume (i) that the propagation of SD can be 
explained even for K — 0, but (ii) that only for K + and appro- 
priate choices of A\j and F the observed two-dimensional SD 
patterns are described correctly, because augmented transmis- 
sion determines the specific spatial evolution of an SD wave on 
the two-dimensional cortical surface. 
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We pose one further constraint on augmented transmission, 
in addition to |A,j| < 1 and \K\ < 1, namely on the choices of 
F. The operators F should be, in contrast to volume transmis- 
sion, of nonlocal nature. These two constraints, summarized as 
weak nonlocal coupling, are the key to incorporating just the 
right amount of flexibility in this framework. The nonlocality, 
on one hand, addresses wiring transmission by describing the 
effect of the structural and functional cortical connectivity on 
SD, and, on the other hand, includes changes in cerebral blood 
flow in response to SD and how these changes are fed back to 
change subsequent neural activity, a feedback process termed 
neurovascular coupling. 

Wiring transmission and neurovascular coupling can include 
nonlocalities in both space and time. For wiring transmission, 
nonlocalities in space are due to long-range lateral cortical con- 
nections, and nonlocalities in time are due to increased open 
probabilities of metabotropic receptors in the range of seconds 
after their activation. The traveling speed of action potentials 
along lateral cortical connections and the ionotropic synaptic 
transmission are, compared to the speed of SD, quasi instan- 
taneous. Nonlocality in space for neurovascular coupling, that 
is, how localized the vascular response is, also referred to as the 
vascular point spread function, depends on which component of 
the vascular system is considered. While changes to capillaries, 
i.e., the smallest component controlled by pericytes 13611 . can 
be considered as a localized transmission in space, compared to 
the SD wave size, changes to arterioles and arteries have wider 
point spread functions. Latencies in the neurovascular coupling 
in the order of seconds are typical and cause nonlocalities in 
time. Moreover, wiring transmission and neurovascular cou- 
pling can act in combination to constitute a nonlocal transmis- 
sion scheme. For example, the SD wave will activate quasi in- 
stantaneously a wide spread increase in neural activity through 
feed-forward and feedback cortical circuitry that induces there- 
fore a spatially global but time-delayed neurovascular response. 

In Sec. 14.51 we introduce a first special case derived from 
Eq. ([T]). This is based on the Hodgkin-Grafstein model of SD 
extended by inhibitor dynamics. This results in a generic model 
with only reaction-diffusion terms (i)-(ii). In Sec. 14.61 we con- 
sider in this generic system emerging patterns and transitions 
between them, in particular, transitions from spiral waves to re- 
tracting wave segments. In Sec. [5] we introduce further cases 
derived from Eq. (HJ including a nonlocal term (iii). 

4.5. Extended Hodgkin-Grafstein model 

The first reaction-diffusion model of SD, suggested by Hodg- 
kin and Grafstein 12711 . was a single species model. It is a spe- 
cial case model that can be derived from our framework for 
e — 0. Then it follows that rate of change d,v = and thus the 
inhibitor v is simply a constant parameter. The activator agent 
was suggested to be extracellular potassium We continue 

to use u as the activator. Moreover, we want to note that this ac- 
tivator can also be a combination of quantities having a positive 
feedback loop, as shown in Fig. [4] without restricting the main 
idea of the Hodgkin-Grafstein model. The only assumption is 



that the activator rate equation /(«, v) is a cubic function in u 



du 
~dt 



(6) 



This model describes a bistable state in cortical ion homeostases. 
The smallest and the largest roots of the cubic function are sta- 
ble fixed points, the lower one being identified as the physiolog- 
ical state, the other one as the pathological state. Furthermore, 
the middle root is an unstable fixed point, i. e., a threshold. The 
Hogkin-Grafstein model equation (|6]l can be solved analytically 
and this provides detailed mathematical insight in the mech- 
anism of how the pathological state invades the physiological 
state (see for example the textbook 13711 ). 

We can choose suitable inhibitor dynamics from the Fitz- 
Hugh-Nagumo equations ll38l . l39ll to extend the Hodgkin-Graf- 
stein model, i. e., to obtain a model with s + 0. This way we 
couple Eq. (O to another ordinary differential equation 



i dv 



— u + (3 — yv, 



(7) 



where /3 and y are two parameters. These parameter determine 
the threshold, which in the single species model (Eq. ©) was 
determined by v. Usually, y is set constant and then /? is solely 
identified as a threshold parameter. 

The FitzHugh-Nagumo equations d6]l-(|7]i are a paradigmatic 
model for excitable media. An excitable medium has the ca- 
pacity to propagate a sustained wave as u and v concentration 
profiles with a defined wave front and trailing edge through 
a spatially extended system. Whereas the Hodgkin-Grafstein 
model can just account for front propagation but not for recov- 
ery defining the dynamics in the trailing edge. 

The biophysical interpretation of the variables and terms in 
the FitzHugh-Nagumo equations ©-(Hi is originally based on 



membrane dynamics in a single neuron 138113911 . because these 
equations were derived from the Hodgkin-Huxley model of ac- 
tion potentials. The interpretation of Eqs. ©-(Ell as a model 
for SD is thus very different. In order not to cause confu- 
sion, we call Eqs. ©-(|7]i an extended Hodgkin-Grafstein model 
in this study, although a great body of knowledge about this 
generic system is published under the name FitzHugh-Nagumo 
model (or also under the name Bonhoeffer-van der Pol model 
1 3 EH 01). 

The interpretation of the terms on the right hand sides of 
Eqs. ©-dTJi can be directly read from Fig. [4] In Eq. (O the first 
"+m" on the right hand side provides the positive self-feedback 
loop (dashed and dotted loops), whereas the nonlinear cubic 
term restricts this autocatalytic process and provides a thresh- 
old. The "— v" comes from the inhibitory effect of the ion pumps 
(arrow with circular head between boxes marked as ion pump 
and ion concentration). The main other important terms (as y 
can be set to a constant in the range < y < 1 without changing 
the generic behavior) are the "+m" and the (3 on the right hand 
side of Eq. (0. The " + u" describes positive pathway to the 
ion pump activity when the ion concentration is out of balance 
(solid arrow with pointed head starting from box marked ion 
gradient and terminating at box ion pump), and the parameter 
sets, as already mentioned, the threshold. 
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4.6. Excitability probed by dynamics of functional block 

The extended Hodgkin-Grafstein model of SD as a generic 
reaction-diffusion model of activator-inhibitor type can explain 
the emergence and transition between different spatio-temporal 
wave patterns in two spatial dimensions, for example from re- 
entrant SD waves revolving around a functional block to retract- 
ing waves segments, by changing parameters in Eqs. ©-(|7]i. 
This in turn can be understood as a translation in the parame- 
ter space of this model. If we set y to a constant value in the 
range < y < 1, there remain only two main parameters in 
this set of equations, namely the time scale ratio s of inhibitor 
and activator dynamics and the threshold [}. Both parameters 
are closely tied to the excitable behavior, that emerges in active 
media. Therefore, the parameter plane (s,/3) can be viewed as 
the excitability plane, see Fig. [5] in which certain bifurcation 
lines mark the emergence of specific spatio-temporal patterns. 

One axis of the excitability plane is spanned by the param- 
eter e. This parameter appears in Eq. (|7]). It is also explicitly 
given in the model equation (Q]), indicating that s is indepen- 
dent on the specific model derived from the framework. It is 
usually a small positive number (0<e«;l), which means that 
the recovery processes (arrows with circular head in Fig. |4j) are 
slower than the autocatalytic excitation caused by positive feed- 
back loops (dashed and dotted arrows in Fig. In Fig. [5] the 
horizontal axis of the excitability plane is spanned by the neg- 
ative logarithm of the time scale ratio e, called time scale ratio 
index - lne. To give an example, if s = 0.1, the kinetic rates of 
the inhibitory processes are by one order of magnitude slower 
than those of the activator and the time scale ratio index is ac- 
cordingly 1 . The other parameter, that spans the vertical axis 
of the excitability plane in Fig. [5] is the excitation threshold, 
which depends on the specific model derived from the frame- 
work. In the extended Hodgkin-Grafstein model /3 can be cho- 
sen as a threshold parameter, if y is set constant, for example 
for Fig.[5]we choose y = 0.5. The units of /3 are rather arbitrary, 
but a characteristic value is y3 = 0.6 (if y = 0.5) at which the 
activator-inhibitor dynamics becomes oscillatory for small s. 

In the parameter plane in Fig. [5] certain bifurcation lines 
mark the parameter values at which specific patterns emerge. 
These pattern are also seen in the experimental SD systems by 
monitoring the functional block (spiral tip dynamics). In fact, 
the dynamics of the functional block can be used to probe the 
tissue excitability and locate the occurrence of these patterns 
in the excitability plane in Fig. [5] For example, it was shown 
that under standard experimental condition (Sec. 12. II ) the func- 
tional block of retinal SD waves is an one-dimensional curved 
line (inset (a) Fig. [5]), which itself performs a complex rota- 
tion [22]. This is caused by a combination of effects. Essen- 
tially, a long trailing edge (indicating small e) causes a large 
refractory zone that the spiral tip, which has a high tendency 
to curl in (indicating small threshold p), must avoid leading 
to complex re-entrant manoeuvres. Similar dynamics of linear 
functional blocks is obtained in computer simulations of the ex- 
tended Hodgkin-Grafstein model (i. e., FitzHugh-Nagumo equa- 
tions ©-(|7|) in a certain parameter regime close to the black 
full circle in Fig. [5] 1 43] . This regime is bounded from one side 
by a bifurcation line dC (dotted line) beyond which complex 



time scale ratio index -In e 
1 2 3 




Figure 5: Parameter space of an generic active medium (FitzHugh-Nagumo 
equations adopted from (43ll ). with parameter regions separated by bifurcation 
lines dC, dM. dR^, and dP that mark sudden changes in the spatio-temporal 
patterns from complex, meandering, rigid rotating spiral patterns, to the left of 
dC, dM and dRa,, respectively, to the propagation boundary BP. To the left 
of dP the medium is non-excitable. The inset (a) is the predicted location of 
retinal spiral-shaped SD waves with a linear functional block (solid line) in the 
regime of complex re-entrant patterns. 



rotation of largely linear functional blocks is found (to the right 
of dC). On the other side of dC, simpler interactions between 
spiral tip and refractory zone occur and the functional block is 
described by a meandering pattern. This region is bounded by 
dM (dashed line). Between yet another bifurcation line, dR<x, 
(solid), and dM occur reentrant patterns with a perfect circu- 
lar functional block (also called rigid rotation), indicating only 
weak or no interactions of the re-entrant wave with its refrac- 
tory trailing edge. At least, in the narrow zone dR m and dP 
(shaded gray) only refracting wave segments occur. 

5. Control and modulation by weak nonlocal coupling 

Connol and modulation of reaction-diffusion patterns by 
weak nonlocal coupling are the key features of the framework 
that provide enough flexibility to explain the experimentally ob- 
served data. In this section, we sketch principal mechanisms 
of nonlocal transmission that lead to the driving force and are 
hitherto hidden behind the nonlocal operator described by the 
terms AjjF. This augmented transmission can explain two- 
dimensional wave patterns observed experimentally. The emer- 
gence of these patterns depends on this augmented transmis- 
sion capabilities and would in this form not occur in a pure 
reaction-diffusion system of activator-inhibitor type. Many of 
these principles have been investigated in other active media in 
detail in previous studies both experimentally and theoretically. 
For example, two-dimensional pattern formation in the chem- 
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ical Belousov-Zhabotinskii reaction was studied with external 
feedback loops. Therefore we only outline these principles and 
refer to the corresponding studies for details. We begin with de- 
scribing changes in excitability as translation in the excitability 
plane, as shown in Fig. |5j along which the excitable state can 
cross bifurcation lines leading to an abrupt change in the two- 
dimensional wave pattern (Sec. 15.11 ). Then, we consider the 
stabilization of wave segments to account for sustained prop- 
agation (Sec. 15.21 ). and, finally investigate re-entrant wave pat- 
terns in this framework with further reference to clinical data of 
persisting migraine aura without infarction (Sec. 15. 3b . 

5.1. Control by long-range connections and time-delayed feed- 
back 

In this subsection, we describe the effect of long-range con- 
nections and time-delayed feedback on reaction-diffusion sys- 
tem of activator-inhibitor type. This was reported in detail in 
recent studies 11441 13311 . Therefore we can limit our description 
to issues that are relevant to link this theoretical work to the ex- 
periments of SD waves in retina (Sec. 13. Q , in particular, to the 
transition from spiral-shaped SD waves to SD wave segments 
by administering drugs (see also Sec. 14.6b . In the outset of this 
section, we consider the efficiency of general control schemes 
causing translations within the excitability plane. 

Let us assume, for the sake of argument, that we can delib- 
erately control the values of the parameters s and f3 that span 
the excitability plane shown in Fig. [5] Furthermore, the system 
is in the state marked by the black full circle in Fig. [5] and the 
control goal is to reach a target state being weakly susceptible 
to SD. In this target state, all SD waves with open ends shall 
retract or even collapse so that SD occurs only transiently. The 
whole regime of states being weakly susceptible to SD is lo- 
cated in the excitability plane to the left of dRoo- Between dRco 
and dP only retracting wave segments can occur (but see also 
Fig. [7] because the location of the rotor boundary actually de- 
pends also on the size of the wave segment), and near but left to 
dP in Fig.[5]only collapsing waves occur. Then the seemingly 
simplest way to reduce tissue excitability is moving along ar- 
row (A) in Fig. [5] The system then moves away from a location 
in the parameter regime of complex re-entrant patterns revolv- 
ing around a linear functional block (Fig. [5] inset), i.e., from a 
location to the right of dC (dotted line), towards a state in the 
regime being only weakly susceptible to SD. 

The efficiency of any control method that moves the sys- 
tem within the excitability plane, for example along arrow (A) 
(Fig. |5]), should be defined by some optimization criterion. Path 
(A) is clearly better than the path along the direction of the ar- 
row (B), because the latter path will not reach the target state, 
although the excitability might seem to decrease moving along 
the direction of path (B), if the dynamics of the functional core 
is used to probe excitability. The core motion becomes sim- 
pler to the left of dC (dotted) and again simpler crossing dM 
(dashed): from complex to meandering to rigidly rotating spiral- 
shaped SD waves, respectively (Sec. 14.61 ). This sequence of pat- 
terns can erroneously indicate that eventually the system might 
become weakly excitable as the same bifurcation lines need to 
be crossed before the target state can be reached. This generic 



behavior of complex systems poses a problem, if the monitored 
functional core dynamics are used to probe the effect of anti- 
migraine drugs. The same holds if even simpler wave param- 
eters such as the velocity are used to probe the effect of anti- 
migraine drugs 1451 146114711 . Moreover, since a parameter plane 
has no metrical structure, other measures than comparing paths 
that reach the target state to those that do not are not available, 
in particular there is no privileged shortest path. Without fur- 
ther constraints there exists an infinite number of possible target 
states in the weakly susceptible regime. 

Let us now assume that excitability can only be controlled 
via nonlocal transmission described by the operator F, in par- 
ticular, via the influence of long-range connections or time- 
delayed feedback. The control parameter is then the coupling 
strength K in Eq. [3] We limit the nonlocal transmission types 
to the simplest coupling schemes with only one element of the 
coupling matrix Ay being one and the others zero. This leads 
to four coupling schemes, two self-coupling schemes with i = j 
and two cross-coupling schemes with ; + j. Other coupling 
schemes, for example diagonal coupling with A, ; = 1 for i — j 
and zero otherwise, are also possible and can be treated in anal- 
ogy, but the limitation to only one non-vanishing component A, ; 
provides the four simplest algebraic constraints that allows us 
to evaluate target states in the weakly susceptible regime of the 
excitability plane (Fig. [5]), if the system is outside this regime. 

Let the matrix element not equal to zero be denoted as A ws 
with w, s e {u, v}, so that the only augmented transmission is the 
following. The state from species s is used to create a nonlocal 
driving force F[s] that is fed back to species w. 

We first consider the effect of long-range connections on 
a planar SD wave extending in the y-direction, so that only 
the propagation along the ^-direction needs to be considered. 
Moreover, we introduce the long-range connections only in one 
direction pointing against the direction of wave propagation. 
Without loss of generality, we let the wave propagate in the 
negative x-direction, then this long-range connection leads to 
the additional term in the rate equation of specie w 



F[s(x, t)] = s(x -6,t)- s(x, t). 



(8) 



where 6 is the distance of the spatial long-range connection. 
This directed connectivity corresponds to anisotropic nonlocal 
coupling in a two dimensional active medium. The functional 
and structural connectivity of the cortex is realistically mod- 
eled as an anisotropic (and also inhomogeneous) medium due 
to the patchy nature of nonlocal horizontal cortical connections, 
however our motivation to introduce Eq. © was different. One 
reason is to obtain a better understanding of the results obtained 
in previous work lO . where isotropic nonlocal coupling have 
also been considered. Moreover, there is a direct analogy to lo- 
cal time-delayed feedback. In case of local time-delayed feed- 
back the control force F is given by 



F[s(x, t)] = s(x, t - t) - s(x, t), 



(9) 



where t is the delay time. 

The common feature of the nonlocal transmission types in 
Eqs. ©-(O is that the front of the reaction-diffusion wave de- 
scribed by Eqs. ©-(|7]) receives an additional driving force from 
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the cortical steady state that lies ahead of the wave (in space or 
time). The main idea, first introduced in Ref. 13311 . is to replace 
the nonlocal terms, that is, the time-delayed term in Eq. (0, 
and the space shifted term in Eq. ©, by the fixed point values 
s(x, t — r) — s(x + 6, t) = s* of species s for this is the state ahead 
of the propagating wave. For example, for the coupling scheme 
A vu this special case derived from Eq. (Q]i with Eqs. ©-(|7]i be- 
comes 



Migraine SD patterns 
(b) 



du 



v + V 2 u 



s — = u + p — yv + K(u — u), 
at 

which can be transformed by simple algebraic operations 

e = (l-K)s 

p = (fi + Ku*)/0.-K) 

y = y/(l - K) 



(10) 
(ID 



(12) 
(13) 
(14) 



to the reaction-diffusion system of activator-inhibitor type Eqs. 
©-(O with effective parameters s, {3, and y. 

The effective parameters s(K), /3(K), and y(K) in Eqs. (TT~2T >- 
( fT~4-b describe a motion along a defined pathway parametrized by 
the coupling strength K in the parameter space of the extended 
Hodgkin-Grafstein model. This pathway will cross dRoo at a 
certain point, while another coupling scheme A ws can lead to a 
different pathway that crosses dRoo at a different point for a dif- 
ferent value of the coupling strength K. In fact, it can be shown 
that the cross coupling scheme A uv leads also to Eqs. (TT~2T> - (TT~4T> , 
only with a change in the sign of K for the latter. Therefore, the 
cross coupling schemes A vu and A uv do not differ with respect 
to their effect. Note, however, that whether they differ with re- 
spect to efficiency must be decided in a dimensionalized model 
and there the ratio (sA vu )/A uv must be used to compare the rela- 
tive coupling strengths of the cross coupling schemes, because 
the physical units must be the same (in analogy, the ratio of the 
diffusion coefficients D v /D u = sd, see Sec. 14.3b . 

By changing the coupling strength K all three parameters 
are varied, therefore we cannot directly illustrate the shift in the 
excitability plane (s,/3) in Fig. [5] which is a cross-section of the 
full parameter space taken at y — 0.5. How the effective pa- 
rameters reduce susceptibility to wave propagations described 
by either Eqs. ([T2l-([T4Ti or equations corresponding to the other 
coupling schemes A uu , A uv , and A vv , was previously investigated 
in detail JH. 

The extended Hodgkin-Grafstein model is based on the Fitz- 
Hugh-Nagumo equation, so this approach describes generic re- 
sponses of reaction-diffusion systems in the form of activator- 
inhibitor type to control by nonlocal coupling types given in 
Eqs. (O-®. We want to stress a different perspective, namely 
that of a control failure and the subsequent change of the cor- 
tical state that was before the failure in a state nonsusceptible 
to SD. An internal nonlocal cortical control mechanism in the 
form of AijF that renders the cortical state nonsusceptible to SD 
and that under pathological conditions is attenuated (K — > 0) 
can drive the cortex into a state where retracting SD wave seg- 
ments (weakly susceptible) occur or even re-entrant SD waves. 





temporarily 
affected area 




Figure 6: Paradigms of SD propagation in migraine. Current view adopted from 
0]: (a) SD wave engulfing a whole hemisphere, (b) A full-scale attack affects 
most of the cortical surface (gray area). SD propagation is assumed to stop at 
the central and lateral sulcus (CS and LS, dotted line), omitting the frontal lobe 
(not shaded). Proposed change in paradigm: (c) localized SD wave segment 
affecting in a full-scale attack only a confined region of the cortex (b). 



This provides a putative common mechanism of SD in migraine 
and stroke iHl whereas the cause of the changed control path- 
way can be very different for the two conditions (cf. Fig. [7] 
insets (a) and (b)). 

5.2. Stabilization of wave segments 

Fig- 0(a) illustrates an SD wave that was triggered at some 
previous time (about lOmzn earlier) in the rearmost portion of 
the human cortex (occipital lobe). Such figures became paradig- 
matic illustrations of SD wave propagation in migraine with 
aura, found in modern textbooks of headache 14811 . This illus- 
tration made its first appearance in a seminal paper spearhead- 
ing the SD theory of migraine more than 20 years ago 0]. It 
suggests that within a single migraine aura attack lasting about 
20min, half the cortex is transiently invaded. Fig. [6](b) shows 
the spatial extent of a "full-scale attack" Jj]. The wave stops in 
the middle, at the central and lateral sulcus (CS and LS, dotted 
line), and omits the frontal lobe in these illustrations. The sud- 
den stop remains unexplained, although there are speculations 
in the literature suggesting "some striking metabolic difference 
of the two regions" 114911 . 

The term full-scale attack refers to the maximal spatial ex- 
tension of an SD wave in migraine. In many attacks, the SD 
wave will stop before. Even taking into account less severe at- 
tacks, in which SD stops earlier, it has been suggested that this 
paradigm illustrated in Fig.|6](a)-(b) must be modified, namely 
that the maximal extent of invaded cortical tissue in a full-scale 
attack is far more limited (Figj6| (d)) caused by a localized 
SD wave segment (Fig.[6](c)) 12 ill . The instantaneously acti- 
vated area by SD in the human cortex is thus confined in its 
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width. From a pattern formation point of view taken within the 
framework of Eq. (d}-©, this is quite a different picture reveal- 
ing the presence of an additional nonlocal feedback mechanism 
{K + 0). In the next paragraph we put forward arguments in 
favour of the pattern in Fig. |6](c)-(d). The last of these argu- 
ments hints at the nature of a possibly involved nonlocal feed- 
back represented by the operator F. 

If we assume that an SD wave engulfs on its course half the 
cortex as shown in Fig.|6](b), then only a subset of this activa- 
tion results in sensory, e.g., visual, awareness This seems 
rather unlikely as the cellular depolarisation during SD is the 
most dramatic cortical event known. Therefore, the alternative 
is more likely, namely that SD in migraine occurs as a localized 
wave segment that within a full-scale attack affects cortical ar- 
eas along certain, though variable, paths (Fig. [6] (d)). This is 
supported by a study investigating the sequence of neurologi- 
cal symptoms in migraine 15111 . Moreover, similar propagation 
patterns were observed for the gyrencephalic feline cortex |52fl : 
while the primary SD wave engulfed one cortical hemisphere, 
corresponding to Fig. [6] (b), subsequent secondary waves re- 
mained within the originating gyrus, similar to Fig.|6](d) though 
more fragmented. Finally, there is one study using functional 
magnetic resonance imaging (fMRI) in humans during a mi- 
graine attack 01 ■ Although the spreading activity pattern ob- 
served by fMRI seems to support the current paradigmatic view 
(Fig. 0(a)), the data is still inconclusive, as noted by Wilkinson 
[50], because most of the observed spreading activation could 
represent network activity through feed-forward and feedback 
circuitry but not a depolarization wave. If this is correct, the 
global activity in the neural circuitry, while not being part of 
SD, can modulate cortical susceptibility to SD. This kind of 
global activity is, as we will show, exactly what is needed to 
understand stable propagation of wave segments 

Suggesting that an SD wave in migraine takes the shape 
of a wave segment leads to one main problem. Wave seg- 
ments are usually not solutions to reaction-diffusion systems 
of activator-inhibitor type, and therefore such forms, chosen 
as an initial condition, transform quickly. An open end of an 
SD wave segment either curls in to form a spiral-shape wave 
(Fig. [5] inset). This usually happens if cortical susceptibility 
to SD is rather high. To be precise, open ends of sufficiently 
large wave segments curl in, if the parameters s and /3 of the 
system lie towards the side of the rotor boundary dRoo where 
the value of the threshold (3 is lower and the value of the re- 
covery rate is slower (larger time scale ratio index - In e) than 
it is on the nearest point lying on dRco- On the other side of 
the boundary dRco, wave segments independent of their size re- 
tract quickly. The observed retinal SD wave patterns (Fig. |2]i 
retract and vanished within a few millimeters (or, in temporal 
units, a few minutes). This is in agreement with wave segments 
that are based on a reaction-diffusion mechanism as described 
by Eqs. ©-(Ell and simulated in the weakly susceptible regime. 
They retract within the same order of magnitude (note that for 
this estimation, Eqs. ©-(Ell need to be dimensionalized, see 
Ref. ||44fl.) Similar results were obtained in a descriptive math- 
ematical model considering the motion of curves with free ends 
1 2411 . In this geometric model, spatially confined SD wave seg- 
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Figure 7: Scheme of the classification of excitability according to spatio- 
temporal patterns (re-entrant/spiral, particle-like, retracting, collapsing) in 2D 
along the dashed-dotted line in Fig. [5] The insets show wave front at two in- 
stances of time (fi : gray, ti'. black) and trajectory of open wave ends (dashed), 
(a) The scheme illustrates a shift in the threshold by an attenuated (K — > 0) 
long-range coupling or time-delayed feedback according to Eq. 1131 (b) The 
scheme illustrates a shear of the control line from dash-dotted to solid as global 
feedback (Eq. ll7t is attenuated (K — > 0). At K = K m (solid line) a saddle-node 
bifurcation occurs. 



ments were first suggested to cause aura symptoms, and it was 
estimated that a full-scale migraine aura attack hardly spreads 
more than several centimeters, and this maximal distance is 
only obtained if parameters are correctly adjusted. 

Adhering strictly to the definition of the rotor boundary as a 
boundary for spiral waves, this boundary depends on the wave 
segment size S and is called dR. At any point to the right of 
dRoo in Fig.[5]only SD wave segments above a critical size curl 
in at their two open ends to form a counter-rotating double spi- 
ral wave, while those below this size retract and the wave van- 
ishes. This can readily be deduced from the fact that the two 
attracting final states coexist in this parameter regime, namely 
the one with a counter-rotating double spiral wave and the ho- 
mogeneous steady state. There is a separatrix manifold in the 
phase space separating the two basins of these two attractors 
and thereby defining a perturbation threshold. Usually, on this 
separatrix manifold a fixed point exists that is stable with re- 
spect to the flow parallel to this manifold. This fixed point is 
a saddle point in the full phase space and corresponds to an 
unstable solution, a so-called particle-like wave of finite size 
revealing a naturally shaped wave segment that only retracts or 
grows when perturbed. 11541 15511 . 

The wave size S can be defined as the surface area of the 
wave segment where a certain elevated activator concentration 
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u e is exceeded 



S(t) 



J ©(« 



(r, t) - u e )dr, 



(15) 



where is the Heaviside function and r denotes the location on 
the two-dimensional cortical surface. Since the wave size S is 
a derived system quantity, it is not represented in the parameter 
plane in Fig. [5] and therefore dR cannot be represented either. 
Instead, the boundary dR,x, is marked as the limit of infinitely 
large wave s egmen ts that neither retract nor grow, so-called crit- 
ical fingers |5J 



5711 . At dReo, the cortical susceptibility to spiral 



SD waves approaches its lower bound, because arbitrarily large 
segments retract below this boundary and therefore waves with 
open ends can not re-enter and will eventually disappear. 

It was shown in a chemical model system of reaction-diffusion 
waves experimentally and with a mathematical reaction-diffusion 
model 11541 15511 that particle-like waves can be stabilized if the 
wave size S (t), as given in Eq. ( TT5l ). is used to adjust excitability 
(p by the linear relation 



4>(f) = aS(t) + b, 



(16) 



where a and b are characteristic feedback parameters of this 
glob al feedback control scheme. In this chemical model system 
1 55J the quantity tp corresponds to the light intensity that con- 
trols the excitability of the chemical reactions, and this quantity 
has a similar function as (3 in the FitzHugh-Nagumo equations. 

The global feedback control scheme in Eq. (fTol l can be trans- 
lated into our framework by setting a — K and b — fio - KS o 
with some reference values f3 , So and replacing in Eq. © the 
constant parameter/? with the controlled threshold /3(f) = <p(t) = 
aS (t) + b so that 



/3(f)=/3 + K(S - So). 



(17) 



For K — 0, the uncontrolled FitzHugh-Nagumo inhibitor rate 
equation with /3 — j3o is obtained. The operator F is then 



F[u(x,t)] =S(t)-S 



(18) 



and all components of the matrix fulfill Ay = except for A vu = 
1. So and K are, like a and b in the notation of Eq. (fTol l, the 
two characteristic feedback parameters of this global feedback 
control scheme. In addition to these parameters, the value of /3o 
plays a central role for this global control scheme. 

The meaning of So, K,f3o, and further related quantities can 
be directly read from the plane in Fig. [7j This is an illustrat- 
ing scheme where the vertical axis is spanned by the controlled 
threshold [3 and corresponds in principle to the dash-dotted line 
marking a section in the parameter space at constant s in Fig. 
[5] The horizontal axis is spanned by the size S of the wave seg- 
ment. Note that we have dropped in the notation of S and /3 the 
explicit time dependence, as this is more compact. Throughout 
this section and in Fig. [7] S and /3 should be read as time depen- 
dent quantities. Both quantities are projections form the infinite 
phase space of the reaction-diffusion system, for S this projec- 
tion is defined via Eq. (15[ and for (3 via Eq. ( TPTT i. Since (3 is a 
quantity derived from S , only a one-dimensional section of the 
(S,f3) plane corresponds to a phase space projection of a system 
with parameters So, K, and fio, as described in the following. 



Without the global coupling, i. e., for K — 0, the value of 
the threshold j3 does not depend on S and it is therefore con- 
stant (J3 = Po). In this case, the evolution of wave segments 
represented in the (S,/3) plane is confined to a horizontal line, 
e. g., the thin dotted horizontal line in Fig. |7j As described in 
previous Sec. 15.11 the value of the threshold of the free system 
(K = 0) can be in three characteristic regimes: (i) either it is in 
the regime where all wave profiles collapse independent of size 
and shape (j3o > Pap), (ii) or where all wave segments retract 
(j3op > fio > y6,3R„), (iii) or where only wave segments below a 
critical size, marked by the solid thick black line dR, retract and 
if having a larger size grow to form a spiral pattern (/3q < fiapj). 
Note that a horizontal line located at /3 — @o indicates the ab- 
sence of a global control scheme in Fig. |7] but j3o can still be an 
effective parameter as given in Eq. ( fT3l l. that is, /3o can be, for 
example, being controlled by time-delayed feedback or long- 
range coupling (inset (a)). 

With global coupling, i. e., for K + 0, the threshold de- 
pends on S and the evolution of a wave segment represented 
in the (S,{3) plane is then confined to an inclined line, e. g., the 
thin solid and dash-dotted lines with arrows. With increasing 
coupling strength K, the inclination with respect to the horizon 
increases and the line of controlled wave segment evolution is 
sheared along the vertical line S = So- So there is one in- 
variant pivot point (So,f3o). If this point is located at certain 
positions in the (S,/3) plane, and dR is a convex boundary of the 
re-entrant wave regime, there exists a value K — Ksn at which 
the inclined line becomes tangent to dR. For values K > Ksn 
(dashed-dotted line, inset (b)), all wave segments retract irre- 
spectively of their size, though the retraction will be slower than 
for any uncontrolled system because the threshold is lowered as 
the wave segment size shrinks. At K = Ksn, particle-like wave 
solutions are born in a saddle-node bifurcation. 

To obtain a saddle-node bifurcation at Ksn, the invariant 
pivot point (So,Po) must be located in either of two regimes: in 
the gray shaded area below the value of /3sr x , i. e., So < dP(fio) 
and Po < PdR„ (shown in Fig.|7ji, where dP(J3o) denotes the pro- 
jection of S onto dP along /3o', or in the rectangular area defined 
by So > S c „ and f3o > PdR*,- If the invariant pivot point (So,/3o) 
is located outside these regimes, again two regimes should be 
distinguished. If (So,fio) is in the rectangular area defined by 
So < Sen and (3o > y6a««,, there exists neither a value of K for 
which the inclined line becomes tangent to dR nor will the in- 
clined line of controlled wave segment evolution intersect with 
dR for K > and therefore particle waves cannot be stabilized. 
If (So,/3o) is in the white area below the value of f3gR x , i.e., 
So > dP(J3o) and [3o < y6as„ there is no critical value of Ksn 
above which global control successfully blows up the basin of 
attraction of the homogeneous steady state to invade the whole 
phase space. 

It is noteworthy to mention that this control adds a distinctly 
new character to the spatio-temporal patterns by merely chang- 
ing the stability of a solution that exists in the uncontrolled sys- 
tem. For example, for K slightly smaller than Ks n, as indicated 
by the solid line in Fig. |7j an initial wave segment above a crit- 
ical wave size S cr u will grow while the value of the controlled 
threshold j3 will simultaneously be adjusted by the global feed- 



12 



back to the final value /T at which a stable particle-like wave 
solution of size S * is obtained. This particle-like wave is also a 
solution, though unstable, of the uncontrolled system Eqs. ©- 
O with /Jo = f>* ■ If, however, the initial wave segment size has 
a value below the critical wave size S cr u, the wave segment will 
retract and once vanished the controlled threshold f3 reaches the 
value Pi, for the homogeneous state of the cortex. 

There are cases where K is too large and the control be- 
comes too hard. For example, for K = 0.1824, y — 0, s — 0.04, 
and /?/, = /?o - KSq = 1.721, we observed that particle-like 
waves can start to temporally breathe, that is, they change their 
size in an oscillatory fashion and reach the final target size S * 
only after a long transient of several oscillations (> 50). From 
this we infer that the real parts of the eigenvalues of the stable 
focus node, which corresponds to the stable particle-like wave 
solution (full black circle), become close to zero. A Hopf bi- 
furcation that leads to stable breat hing pa tterns, as observed in 



Stroke SD patterns 



semiconductor nanostructures 115 8l 15 9116011 . which are based on 



similar dynamics, have, however, not been found. 

It will be an important future task to relate concepts and 
quantities illustrated in Fig. [7J to SD in human cortex in or- 
der to develop new therapeutic methods that make use of these 
concepts. That Fig. [7] correctly describes the control of corti- 
cal susceptibility to SD in migraine is currently only supported 
by various correct predictions of such a model. Beside those 
mentioned in the outset of this paragraph and those reported in 
previous studies 11241 roll l23l 12111 . we note one additional pre- 
diction that shows the unifying character of this framework. 
So-called clinically silent migraine aura, that is, an aura that 
migraineurs are subjectively unaware of because they do not 
experience any neurological disturbance, might simply be ex- 
plained by a quickly retracting SD wave patterns occurring for 
values K > Ksn (dash-dotted line). Likewise, it is known that 
clinically silent epilepsies are caused by seizure activity that 
does not break away from a focus. Our framework would still 
predict that the control mechanism is switched on to suppress 
propagation, which must result in a blood flow 'fingerprint' of 
SD, which indeed is observed by non-invasive imaging 111 411 . 



5.3. Re-entrant SD: a functional definition of tissue at risk 

The re-entrant SD wave pattern observed in the experimen- 
tal stroke model (Sec. l3.2t is illustrated in Fig.[8](a): an anatom- 
ical block of nonexcitable tissue (black) develops surrounded 
by a ring-shaped zone becoming susceptible to SD (gray). The 
experimental procedure (Sec. 12. 2\ does not necessarily cause ir- 
recoverable damage in the form of an infarct lesion in the core 
region, but there is still an anatomical block for a large central 
area is not affected by the cycling. 

The focus in this section is on the surrounding ring-shaped 
zone (gray, Fig. [8] (a)). Its dynamical features are of crucial 
clinical interest, because the aim of stroke therapy is salvage of 
this tissue (Fig. [8](b)). At the outset of this section, we con- 
sider the case that the anatomical block is an infarct core. In 
fact, when infarct lesions in the cerebral cortex are experimen- 
tally induced by vessel occlusion, re-entrant patterns similar to 
those reported here emerge lfl7ll . By the end of this section, 
we will then reconsider that the anatomical block induced with 




Figure 8: (a) Cortical state during a localized graded reduction in blood flow. 
Below a critical threshold of perfusion, an infarct lesion develops (black) sur- 
rounded by tissue becoming susceptible to SD waves (gray). These waves can 
spread in various patterns, e.g. radial outwards from the infarct core and cycling 
around it. (b) The infarct tissue may eventually include much of the surround- 
ing tissue or be limited to the initial infarct core inset. 



our experimental procedure is likely to be only nonsusceptible 
to SD but also salvageable, in particular, we compare the ob- 
served pattern with strikingly similar observations in persistent 
migraine aura without infarction. 

The cortical mechanisms by which the ring-shaped zone 
(gray, Fig. [8] (a)) is formed are not yet fully established and 
consequently there is also not a clear definition of how far this 
zone actually spreads out. We suggest that the mathematical 
framework outlined in Eq. Q]can complement current clinical 
definitions of this spatial structure. The zone surrounding an 
infarct core is either called penumbra j62ll . denoting a concept 
based on a decreased blood flow range, or an alternative con- 
cept is based on the tissue at risk of infarction being potentially 
salvageable. In stroke outcome, the infarct tissue may eventu- 
ally include much of the surrounding tissue (Fig. [8] (b)) or be 
limited to the initial infarct core (Fig.[8](b) inset). 

The tissue-at-risk (TAR) concept is clinically most relevant, 
but it does not formulate a distinctive definition. This concept 
must be substantiated by what processes underlie this risk so 
that it becomes clear what constitutes high and low-risk. A de- 
creased blood flow range is one possible substantiation, which, 
if the sole complement, would identify the penumbra as TAR. 
The penumbra is anatomically defined. By contrast TAR can 
also be defined as a functional concept in terms of nonlinear dy- 
namics 12311 . In its original functional form, TAR was suggested 
to extend over the area invaded by transient waves occurring be- 
low the propagation boundary (to the right of dP (Fig. [5]). This 
functional definition did not incorporate a gradient in suscepti- 
bility caused by anatomical heterogeneities. We suggest in the 
following to combine these anatomical and functional concepts, 
which would then attach great clinical importance to re-entrant 
pattern. 

The defining property of re-entrant pattern with an anatom- 
ical block is the heterogeneity. In stroke, this heterogeneity is 
evolving and occurs actually over several scales. On the largest 
scale is the infarct core and the surrounding region. In the in- 
farct core blood supply is below the threshold of energy failure, 
outside it is constrained but the energy state is preserved. On a 
smaller scale occurs a heterogeneous microcirculation 
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A mathematical model that does not incorporate a complete 
characterization of this heterogeneous microcirculation can use 
nonlocal terms to describe the macroscopic spatial connectivity 
in form of a network of local excitable elements described by 
Eqs. (HrQi). In such a model 1 64116511 . it was shown that spatial 
diversity induces wave patterns. Diversity is able to induce a 
transition to an excitable behavior of the net. These patterns are 
most coherent for an intermediate variability strength, an effect 
that is described as similar to stochastic resonance generated by 
additive noise. 

The spatio-temporal pattern reported in Sec. 13.21 are strik- 
ingly reminiscent of patterns reported by migraineurs suffering 
from a rare but well documented subform of migraine classi- 
fied as persistent migraine aura without infarction in the second 
edition of the international headache classification (ICHD-II 
code 1.5.3). In this subform, migraine aura phenomena, which 
are otherwise transient neurological symptoms lasting less then 
30m;«, persist for more than 1 week, maybe months or years. 
To diagnose this subform, there must not be any evidence by 
noninvasive imaging of an infarct and these symptoms must not 
be attributed to another disorder. In Fig. 0(a) a persistent cen- 
tral visual disturbance is illustrated. The scotoma is filled by ge- 
ometric pseudohallucinations of lattice form. While this pattern 
is seen always, even with eyes closed during an acute migraine 
attack the size of central scotoma can decrease and become up 
to 4 times larger (Fig. 0(b)), or a further disturbance can prop- 
agate in one visual hemifield (c) or in both (d). 

In analogy to the stabilization of particle-like waves de- 
scribed in Sec. 15. 21 these stationary patterns might be explained 
by the stabilization of a critical nucleus with global feedback. 
The critical nucleus refers to a spatially confined structure ly- 
ing in the phase space on the separatrix manifold separating 
the basins of attraction of a ring-shaped wave and the homo- 
geneous steady state. It defines a stationary radial symmetric 
perturbation threshold, while the particle-like wave is a travel- 
ing asymmetric perturbation threshold. The stabilization of the 
critical nucleus can probably not be achieved by global feed- 
back as introduced in Eq. ( fT8l l. because the coupling strength 
must be very large. A large inclination of the control line in 
Fig. [7] can be predicted if assuming that the particle-like wave 
size approaches the size of the critical nucleus Sen f° r small 
thresholds. However, if the operator F has a nonlinear depen- 
dence on S , or additional long-range, time-delayed, and other 
nonlocal augmented transmission schemes, including spatial di- 
versity, are considered in combination with reaction-diffusion, 
a rich repertoire opens up and it will be an important task for 
future research to investigate this in the context of migraine and 
stroke patterns. 

6. Discussion and outlook 

We have presented a mathematical framework for cortical 
spreading depolarizations that unifies existing activator-inhibitor 
models as special cases Within this frame- 

work the reaction-diffusion mechanism of SD is extended by 
nonlocal coupling as an integral part of the phenomenon neces- 
sary to describe the emergence and transitions of two-dimensional 
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Figure 9: (a) Persistent central scotoma seen with open eyes. The area of the 
central scotoma is filled by geometric hallucinations of lattice form dimension. 
The pattern in this area remains when eyes are closed. During a migraine attack 
starting from the central scotoma various patterns can occur, (b) Expanding 
the initial affected area in the visual field, (c) waves propagation in one visual 
hemifield, and (d) waves propagating in both visual hemifields. 



SD wave patterns. As a framework it provides a flexible proce- 
dure for analysing SD patterns observed in both experimental 
migraine and stroke models, and clinical data J5, IT]. 

This framework is aimed at providing basic insights needed 
to understand cortical susceptibility to transient SD wave prop- 
agation in terms of nonlinear science. In particular, we focus 
on bifurcations in the dynamics of cortical homeostasis caused 
by cortical circuits and neurovascular coupling represented as 
nonlocal transmission schemes. A long term goal is to pro- 
vide control strategies based on bifurcation theory for both tra- 
ditional pharmacological treatments and biomedical engineered 
devices that intelligently target the occurrence of SD waves. 

A mathematical model that involves only two species with 
activator-inhibitor dynamics is undoubtedly a very crude macro- 
scopic level of cortical homeostasis. Activator-inhibitor dy- 
namics is the minimal requirement that leads to a description of 
SD as a traveling wave pattern with a trailing edge. In two spa- 
tial dimensions, these waves take a characteristic form that re- 
veals much about cortical susceptibility to SD. Yet, little insight 
into microscopic mechanisms of SD is offered by such minimal 
dynamic systems. There are much more sophisticated mathe- 
matical descriptions of SD on a detailed microscopic level that 
try to resolve the mechanism of SD |66[ 67]. In our framework, 
these microscopic mechanism are hidden inside the box labeled 
as "firing rate" (Fig. [4]). 

Proposing a conceptual framework at this macroscopic level 
must be justified by showing that this approach avoids deficien- 
cies that would arise if SD wave propagation is modeled on a 
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microscopic level. From the clinical perspective, the macro- 
scopic scale, in particular, the two-dimensional pattern on the 
cortical surface, is highly relevant. In migraine, the whole spec- 
trum of aura subtypes could depend critically on whether SD 
can break away from a focus point, and, if it does, whether SD 
waves take the shape of a particle-like wave and therefore limit 
the extent of a full-scale attack to certain cortical regions (Fig. [6] 
(d)), or, in the worst case, SD becomes a persistent activity pat- 
tern (Fig. [9]). In stroke, outcome depends on SD as these waves 
are suspected to contribute to the loss of potentially salvageable 
tissue, i. e., tissue at risk of infarction. Since re-entrant pattern 
increase the total number of waves, the emergence of these pat- 
terns need to be understood. None of these questions can be 
addressed in a cellular model of SD including only one or sev- 
eral neurons, let alone the semantic objection against applying 
the term spreading to such models because they are not meant 
to model the spread ITill . 

Microscopic models include several ion channel types in a 
single neuron, and its dendritic tree. These models also con- 
sider surrounding compartments, and can therefore also describe 
changes in intra- and extracelluar ion concentrations J66, 671. 
To obtain the continuum limit of such discrete microscopic cel- 
lular models is by no means straightforward. The problem is 
to match the very detailed and accurate knowledge of micro- 
scopic processes involved in the depolarization cycle of SD 
by an equally detailed description of volume transmission and 
wiring transmission needed to provide the spatial coupling for 
transition to the continuum limit. Mismatching these levels of 
description by merely using a diffusion term as in Eq. (f2]i to 
simulate spatial coupling by extracellular potassium would be 
like attaching the wheels of a carriage on a Mercedes. 

Efforts have been made to address this problem by account- 
ing not only for interstitial ionic diffusion but also for ionic 
movement through a neuronal syncytium of cells connected by 
gap junctions and for cells that are allowed to expand in re- 
sponse to osmotic pressure. From these it was concluded that 
cytosolic diffusion via gap junctions and osmotic forces are 
important mechanisms underlying SD |68]. Yet, to compute 
the spatio-temporal development of SD across the cortical sur- 
face on a centimeter scale and investigate the emergence of re- 
tracting, re-entrant, and stationary SD waves, coupling schemes 
other than local schemes must be accounted for. There are im- 
portant nonlocal transmission schemes in the cortex hat provide 
spatial coupling over a long range or lead to time delays. Ex- 
ample are (i) the functional and structural cortical connectivity, 
in the continuum limit described by neural fields J69(], (ii) the 
cortical energy state maintained by a segmental blood perfu- 
sion from the arteries to smaller arterioles and finally branching 
capillaries leading to a complex hemodynamic response that 
is regulated by smooth muscle cells (arterioles) and pericytes 
(capillaries) 13611 . and (iii) the topology of the network of the 
neuronal syncytium that will be compromised by disseminate 
neuronal injury during ischemia 117 Oil . There are more trans- 
mission schemes, but the effect of these (i)-(iii) are the ones 
we have considered in a first approximation by (i) long-range 
connections, (ii) global coupling, and (iii) spatial diversity. 

Future investigations have to relate the quantitative values 



of the parameters in macroscopic model that successfully de- 
scribe the bifurcations observed in the spatio-temporal char- 
acteristics of SD to the cellular level. This will open up this 
approach for the design of both drug treatments 12311 and engi- 
neered devices utilizing transcranial stimulation 171(1 with the 
aim to target SD intelligently. Migraine and stroke are found 
in specific genetic disorders giving clues to genetic factors that 
hint at differences in the cortical network activity. Such hints 
may provide means to bridge the gap between the macroscopic 
and microscopic level of SD. Similar hints of changed network 
activity are coming from psychophysics describing abnormal 
cortical processing in migraine by concepts like hyper- and hy- 
poexcitability, heightened responsiveness, a lack of habituation 
and/or a lack of intra-cortical inhibition [72, 5(1 ■ Our frame- 
work provides means by which such statements on cortical ex- 
citability in migraine can be investigated, in particular how ab- 
normal cortical nonlocal connectivity changes susceptibility to 
SD. Although it seems tempting to suggest that cortical hyper- 
excitability increases susceptibility to SD or even that neurons 
prone to hyperexcitabilty trigger SD, such a simple relation can- 
not be expected, and a detailed bifurcation analysis will be im- 
portant for the thorough understanding of SD as a nonlinear 
pattern formation process. 
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